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The recently discovered superconductor Cui;Bi2Se3[T] is a candidate for three-dimensional time- 
reversal- invariant topological superconductors 12], which are predicted to have robust surface Andreev 
bound states hosting massless Majorana fermions. In this work, we present an analytical and 
numerical study of the surface Andreev bound state wavefunction and dispersion. We find the 
topologically protected Majorana fermions at fc = 0, as well as a new type of surface Andreev bound 
states at finite k. We relate our results to a recent point-contact spectroscopy experiment [3] • 
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The discovery of topological insulators has generated 
much interest in not only understanding their proper- 
ties and potential applications to spintronics and thermo- 
electrics but also searching for related topological phases 
in new directions. A particularly exciting avenue is 
topological superconductivity [21 IH-[T5]. in which uncon- 
ventional pairing symmetries generate nontrivial super- 
conducting gaps in a similar way that spin-orbit cou- 
pling generates inverted band gaps for topological insu- 
lators. The hallmark of a topological superconductor is 
the existence of gapless surface Andreev bound states 
(SABS) hosting itinerant neutral Bogoliubov quasiparti- 
cles, which are the analog of massless Majorana fermions 
in high energy physics. 

There is currently an intensive search for topologi- 
cal superconductors. In particular, a recently discov- 
ered superconductor Cu2:Bi2Se3 with Tc SK has at- 
tracted much attention [1. It was proposed that the 
strong spin-orbit coupled band structure of Cu2;Bi2Se3 
may favor an unconventional odd-parity pairing symme- 
try, which naturally leads to a time-reversal-invariant 
topological superconductor!^. Subsequently, many ex- 
perimental and theoretical efforts [TSHTO] have been made 
towards understanding superconductivity in Cua;Bi2Se3. 
In a very recent point-contact spectroscopy experiment, 
Sasaki et a/. [3] have observed a zero-bias conductance 
peak which is attributed to SABS and seems to signify 
unconventional pairing jl4j. 

Motivated by this finding, in this Letter we study the 
phase diagram of odd-parity topological superconductiv- 
ity and the resulting surface Andreev bound states in 
doped semiconductors with strong spin-orbit coupling, of 
which Cuj^Bi2Se3 is a prime candidate. We start from a 
k ■ p Hamiltonian which captures the essential features of 
its band structure near the Fermi surface. By studying 
the phase diagram of the k ■ p Hamiltonian as a func- 
tion of band gap, pairing potential, and doping, we es- 
tablish three gapped phases: topological superconductor 
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FIG. 1: a) Side view of a semi-infinite crystal of Bi2Se3. The 
two relevant pz orbitals are shown in the zoom-in view of the 
QL unit cell, b) Bulk and surface bands of the tight-binding 
model for Bi2Se3. pi and /X2 denote two chemical potentials 
where the surface states have, respectively, not merged and 
merged into the bulk bands. 



(TSC), topological insulator (TI), and normal band in- 
sulator (BI). We characterize these phases in a unified 
way by introducing a topological invariant — a general- 
ized mirror Chern number. We find that the odd-parity 
topological superconductivity in both doped BI and TI 
gives rise to surface Majorana fermions with a linear dis- 
persion at k = 0, as expected from the bulk-boundary 
correspondence. However, the SABS of the two become 
quite different at large k. In particular, we infer from 
the mirror Chern number that the SABS in a supercon- 
ducting doped TI must become gapless again near the 
Fermi momentum. This results in a new type of zero- 
energy SABS. To support these findings, we construct 
a two-orbital tight-binding model to calculate the SABS 
dispersion numerically. Finally we relate these results to 
the recent experiment . 

Band Structure: We begin by reviewing the crystal 
structure of Bi2Se3 and the k -p Hamiltonian for its band 
structure. Bi2Se3 is a rhombohedral crystal in which each 
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unit cell is a quintuple layer (QL) of the form Se(A)- 
Bi(B)-Se(C)-Bi(A)-Se(B). Here A, B, C denote three tri- 
angular lattices stacked with offset along the (111) direc- 
tion, which we will define as the 0-axis. Undoped Bi2Se3 
is a topological insulator. Cu doping introduces elec- 
tron carriers into the conduction band with a density of 
~ 10^°cm~^, which results in a small 3D Fermi surface 
centered at F. The band structure near F is described by 
a, k ■ p Hamiltonianpi EO] : 



H{k) 



+ V^k^CTy -|- V(Tz{kxSy — kySx)- 



(1) 



Here cr^ = ±1 labels the two Wannier functions which are 
primarily Se and Bi orbitals on the upper and lower 
part of the QL unit cell respectively (see Fig.l), and each 
has a two-fold spin degeneracy labeled by = ±1. As 
explained in Ref.[2], the form of i?(k) can be deduced 
entirely by time reversal and crystal symmetry consider- 
ations. The physical origin of various terms in ^ will 
become clear from an explicit tight-binding Hamiltonian 
later. While the magnitudes of the parameters m, Vz and 
V in -ff (k) have been obtained by fitting the band disper- 
sion with the angle-resolved photoemission spectroscopy 
measurement [T5] . their signs turn out to be particularly 
important in this work. 

Phase Diagram: Now consider the superconducting 
Cua;Bi2Se3 with the spin triplet, orbital singlet pairing 
proposed in Ref . P] . The mean-field Hamiltonian is given 

by 
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(2) 



Here i?BdG is the Bogohubov-de Gennes (BdG) Hamil- 
tonian; Tx^z are Pauli matrices in Nambu space; A is the 
pairing potential; ^ is chemical potential. 

The Hamiltonian ([2| exhibits three topologically dis- 
tinct gapped phases as a function of the band gap, pairing 
potential and doping. At zero doping = 0) and in the 
absence of superconductivity (A = 0), the system is ei- 
ther an ordinary BI or a TI, depending on the sign of 
m. At finite electron doping, the chemical potential lies 
inside the conduction band: /.t > 0. When the odd-parity 
pairing A occurs in such a doped BI or TI, the system be- 
comes a fully gapped TSC. For the sake of our argument, 
it is useful to note that the TSC phase is adiabatically 
connected to the /i = and A > |to| limit. Fig. 2 shows 
the BI, TI and TSC phases in the /i = phase diagram as 
a function of m and A. The topological phase transition 
between TSC and BI/TI occurs at A = ±m. 

To characterize these three phases in a unified way, 
we introduce a common topological invariant — a gener- 
alized "mirror Chern number" Um- Recall that reflection 
M with respect to the yz mirror plane is a symmetry el- 
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FIG. 2; Phase diagram of fully-gapped odd-parity supercon- 
ductivity in doped semiconductors as a function of band gap 
m and pairing potential A, showing three gapped phases: 
band insulator, topological insulator and topological super- 
conductor. They are topologically distinguished by the mirror 
Chern numbers um- 



ement of the Bi2Se3 crystal point group D-^d- Under mir- 
ror reflection, the band structure (!]) is invariant, whereas 
the pairing order parameter in (12) changes sign. (Mirror 
reflection flips Sy and s^.) This sign change can be com- 
pensated by a gauge transformation A — )■ —A, so that 
the BdG Hamiltonian is invariant under a generalized 
mirror operation M — AItz, M — ~iSx- 



HBdci^x 



njy , njz 



BdG 



-kx, ky, kz)M 



(3) 



In particular, since i?BdG in the kx = plane commutes 
with M, a mirror Chern number um can be defined 
in the same way as for an insulator [21]. Specifically, 
h = i?BdG(fca; = 0) is a direct sum of two subsystems 
h± with mirror eigenvalues ±i respectively: h± = P±h, 
where P± = (1 ^ iM)/2. The Chern number n± for h± 
satisfies n^+ri- — required by time reversal symmetry. 
However, the difference defines the mirror Chern number: 
Um = {n+ - n-)/2 = n+. 

Using nM(BI) = as a reference, we obtain the mirror 
Chern number for the TI and TSC by calculating the 
change of um across the phase transition to the BI. Due 
to the double counting of electrons and holes, the mirror 
Chern number of an insulator defined in Nambu space 
is twice the value of that defined previously [H]. As a 
result, a direct transition from TI to BI at A = changes 
"n-M by two. For A ^ 0, this transition is split into two 
transitions with an intermediate TSC phase, so that each 
transition changes um by one. Therefore we have 



nM(TI) =2nM(TSC). 



(4) 



The fact that TI and TSC have mirror Chern numbers 
of the same sign will play a key role in our analysis of 
gapless excitations on their surfaces. 

Surface States: Consider a semi-infinite Bi2Se3 crys- 
tal occupying z < 0, which is naturally cleaved between 
QLs (see Fig.l). The boundary condition corresponding 
to such a termination in k ■ p theory is 



azib{z = 0) = V(2 = 0). 



(5) 
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As explained in Ref[5], this boundary condition reflects 
the vanishing of the wavefunction on the bottom layer 
((Tz ~ —1) at z — 0; whereas the other termination within 
the QL (not yet reported in experiments) corresponds to 
a different boundary condition a^'i) = —ip at z = 0|22j. 

First we study surface states of Bi2Se3 by solving the 
differential equation H{kx,ky,—idz)4' = Eip subject to 
the boundary condition ([s]) . We obtain two sets of exact 
solutions ?/'±(k||,z) with the energy-momentum disper- 
sion E±{k) — ±vk. 



V'±(k||,z) = e^/'(l,0),® (l,±ie^^). 



(6) 



where I — ^Vz/m is the decay length; is the azimuthal 
angle of ky ; the subscripts a and s denote the orbital Cz 
and spin Sz basis. In order for the surface states to exist 
in a TI, we must have decaying solutions in the —z di- 
rection, which implies w^m < 0. The spin polarization of 
tlj± is locked to its momentum, forming a two-dimensional 
Dirac cone. In contrast, a BI has Vzm > for which no 
surface state solutions exist. 

Next we study surface Andreev bound states of 
Cua;Bi2Se3. We start by solving the BdG equation at 
k|| =0. A Kramers pair of zero-energy eigenstates 
ipa=±i^) with mirror eigenvalues AI = i ■ a is expected 
from the topology and symmetry of H-BdcO^w ~ ^"f' 
ipai^) satisfies a reduced 4-component equation: 



[{ma^ - ivzaydz - p)tz + AcTyTj-Ji/j" (z) = 0. 



(7) 



By multiplying both sides of Eq.Q with t^, it becomes 
clear that ^jj^ is an eigenstate of Ty. The corresponding 
eigenvalue is given by sgn(wz) in order to have a decay- 
ing solution. Eq.Q then reduces to a two-component 
equation in orbital space, which has two independent so- 
lutions: 



(8) 



Here kp is the Fermi momentum in the z direction, given 
hy kp = \J [i^ ~ m? /vz\ k is the inverse decay length 
given by K = ^/\vz\ > 0; 6* is an angle defined by 
e'^ = {m + iy/ fj? — m?)/ ^. We now choose a suitable 
linear combination of and to satisfy the boundary 
condition ([s]) and obtain the wavefunction of SABS: 

^'"(z) cx e'^''{s\TL{kpz - 0),sin{kFz))„ ® 

[(1, -a)s ® (1, 0)^ + isgn(wz)(l, a), (0, 1)^] , 

where r = ±1 denote the particle/hole basis. ij^g{z) is 
particle- hole symmetric: Si/j^ (z) = ^^^{z) up to an overall 
phase (S — SyTyK), and therefore represents a Majorana 
fermion. In the limit m = 0, agrees with the result 
from Ref.[2]. 

Away from k — 0, the Kramers doublet -0^ and ijj'^ is 
split by spin-orbit coupling term in ([ij. To lowest order 
in k, the SABS dispersion ea{k) is linear and given by 



^aik) = avk. The velocity v is obtained from first-order 
perturbation theory: 

A'^ + sgn{vz)Am 

V — V (9) 

A2 -t- sgn{vz)Am + ^j? 

Since A ^ |r7i| < /i in weak-coupling superconductors 
(provided m is not vanishingly small), ([9]) simplifies to 

■5 ~ w • sgn(wz)ATO//i'^. (10) 

While the SABS in previous studies of single-band su- 
perconductors can often be obtained by quasi-classical 
methods from the geometry of the Fermi surface alone, 
this is not the case for multi- orbital systems such as 
Cua;Bi2Se3. Here both the wavefunction and dispersion 
of SABS depend on the dimensionless factor m//^, which 
arises from the orbital character of electron wavefunc- 
tions on the Fermi surface. 

To deduce the behavior of SABS dispersion for large 
fc, we make use of the mirror Chcrn number introduced 
earlier. According to the principle of bulk-boundary cor- 
respondence, the sign of mirror Chern number um deter- 
mines the helicity of surface states as follows. nM < im- 
plies that the branch of surface states at k^ = with mir- 
ror eigenvalues =Fi (i.e., Sx — ±1) moves (anti-)clockwise 
with respect to +x axis at the edge of the yz plane, which 
means that the electron's spin and angular momentum 
from circulation at the edge are parallel. Such a helic- 
ity forces the dispersion of the Sx — +i- surface band 
to eventually merge into the E > bulk quasiparticle 
continuum at a large positive k. 

Given this correspondence between the mirror Chern 
number and the helicity, the SABS in the TSC must have 
the same helicity as the TI surface states, because the 
mirror Chern numbers Q of the two phases have the 
same sign. Moreover, the helicity of Bi2Se3 surface states 
is directly given by the sign of Dirac velocity v (since the 
dispersion is monotonic). Comparing the SABS velocity 
V in (10) with v, we find that sgn({;) = sgn(u) for TSC 
in a doped BI {vzm > 0), whereas sgn(w) = — sgn(w) for 
TSC in a doped TI {vzm < 0). The latter case deserves 
further attention. In order to fulfill the helicity require- 
ment, the two branches of SABS with opposite mirror 
eigenvalues must eventually switch places and thus cross 
E = once again at some finite k (or an odd number of 
times). As long as perturbations do not change the sign 
oi V at k = 0, the second crossing is robust due to the 
mirror Chern number argument. Similar reasoning was 
previously applied to surface states in TIs|2 H l23l I24j . 

To demonstrate the SABS explicitly, we construct a 
two-orbital tight-binding model in the rhombohedral lat- 
tice shown in Fig.l and calculate the SABS dispersion 
numerically. The Hamiltonian is defined as follows: 



H = Hq + H12 



(11) 



^0 = E<,j> ^ocLcja and H12 = E<»Gije2> ^iclaCj" + 
S<jei j'G2> ^2 c|(jCj'Q, describes spin-independent nearest 
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FIG. 3: SABS dispersion for the tight-binding model in which 
a.) m — —0.3 < 0, /ii = 0.6 and h) m — —0.3 < 0, ^2 = 1 
correspond to doped TI; c) m — 0.3 > 0, /xi = 0.6 corresponds 
to a doped BI; d) m — 0, /ii — 0.6 corresponds to a doped 
zero-gap semiconductor. The mirror eigenvalues are displayed 
near each branch of SABS, Note that in a doped TI, SABS is 
twisted with a second crossing near Fermi momentum. 

neighbor hopping within the same layer (to), as well as 
between two neighboring layers within a QL (ti) and 
on two adjacent QLs {t2)- i?soc describes the Rashba- 
type spin-orbit coupling associated with nearest-neighbor 
intra-layer hopping, which take opposite signs on the top 
and bottom layer of the unit cell due to opposite local 
electric fields along the z direction: 

<ij>ei <ij>e2 

where a^j — ^eijk(Rj — Rfc) denote the vectors joining 
nearest neighbors within a layer, and Ri,2,3 are the Bra- 
vais lattice vectors. The last term H'12 describes inter- 
layer second nearest neighbor (^3) hopping within a QL: 

^'12 ~ S<<ielje2>> ^S^ia^ja + h.c. 

This tight-binding model is constructed to reproduce 
the k ■ p Hamiltonian (1) of Bi2Se3 in the small k limit 
as follows: m = 3(ti +t2 + ^3), Wz = 3t2C, and v = |Aa^, 
where a = |ajj| and c = ||(Ri + R2 + Rs)!- Moreover, 
since there is hopping beyond two neighboring layers in 
H, the open boundary condition in this lattice model 
correctly reproduces the desired boundary condition ([s]) 
in the continuum theory [57]. It is to be emphasized that 
this tight-binding model does not aim to describe the 
band structure of Bi2Se3 in the entire Brillouin zone|28]. 

The bulk and surface bands of the above tight-binding 
model are displayed in Figure lb. To include supercon- 
ductivity, we add the following odd-parity pairing term 
in the Hamiltonian: 

hmf = h+ J2 t + ^l^t) + 

<iei,ie2> 



Fig. 3 shows the resulting SABS dispersion in both doped 
TI (m < 0, Fig.3a-b) and doped BI (to > 0, Fig.3c). The 
SABS in both cases have a branch of linearly dispers- 
ing Majorana fermion at = 0, which signifies a three- 
dimensional topological superconductor. Moreover, the 
SABS in doped TI has a "twisted" dispersion with a sec- 
ond crossing near Fermi momentum, as anticipated from 
our earlier mirror Chern number argument. 

The tight-binding calculation sheds light on the origin 
of twisted SABS. The presence of the second crossing is 
a remnant of the TI surface states in the normal state. 
This can be easily understood in the case where electron 
surface states remain well-defined at the Fermi energy 
in the normal state, e.g, at the chemical potential fj,i 
in Fig.lc. The surface states at k and — k have oppo- 
site mirror eigenvalues, whereas A only pairs states with 
the same mirror eigenvalues. Due to this symmetry in- 
compatibility, surface states must remain gapless even in 
the presence of such an odd-parity pairing. Therefore, 
the Majorana fermion SABS at fc = smoothly evolves 
into TI surface states near Fermi momentum. Indeed we 
find that the decay length of SABS at and beyond the 
second crossing is comparable to that of the electron sur- 
face state given by |w2|/|to|, which is much shorter than 
the typical decay length of SABS given by jw^l/A. This 
is further confirmed by an analytic calculation using the 
k ■ p BdG Hamiltonian ([2| . We search for an E = solu- 
tion at a nonzero fco and find ko is the nontrivial solution 
of 

2Re{W) + ^{~l + \W\^) = 0, (12) 

where 

^{vkoY + {^ + iEpY' 

For A = and to < 0, the solution is given by fco = m/'^j 
which represents the TI surface state at Fermi energy. 
For A 7^ 0, the solution is perturbed: fcg = ^(1 — ^^) 
to leading order in A. No solution exists for to > 0. 

What happens if there are no surface states at the 
Fermi energy in the normal state of a doped TI? For 
example, at the chemical potential pi2 in Fig.lc, the TI 
surface states have merged into the bulk|5S]. Remark- 
ably, after turning on the odd-parity pairing (Fig. 3b), 
the SABS still has the second crossing, as required by 
the mirror Chern number. The resulting gapless SABS 
near Fermi momentum have significant particle-hole mix- 
ing and therefore cannot be interpreted as unpaired TI 
surface states. This represents a new type of surface An- 
dreev bound state, which originates from the interplay 
between band structure and unconventional supercon- 
ductivity. Such SABS defy a conventional quasi-classical 
description and are worth further study. 

Finally, we relate our findings to the recent point- 
contact spectroscopy experiment on CuxBi2Se3[2j. A 
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zero-bias conductance peak with a width of 0.6mV was 
clearly observed below 1.2K, which indicates the exis- 
tence of SABS. A natural candidate is SABS of the fully- 
gapped odd-parity topological superconductor phase pro- 
posed in Ref.|2]- In this Letter, we found that the SABS 
of TSC in a doped TI have unusual and topologically 
robust features. As shown in Fig. 3a and b, the second 
crossings of SABS at finite k along all in-plane directions 
make the density of states at zero energy finite, in sharp 
contrast to the vanishing density of states at /c = char- 
acteristic of the linear dispersion in 2D as found in a 
doped BI (Fig. 3c). Moreover, for small m//Lt the SABS 
dispersion is quite flat as shown in Fig. 3b and 3d. These 
features may naturally explain the experimental obser- 
vations, and will be studied in a future work. 

While the main focus of this Letter is Cua;Bi2Se3, we 
end by discussing the applicability of our findings to 
a large class of superconducting doped semiconductors 
with inversion symmetry. Potential candidates include 
BiaTealSS] and TlBiTesISI] under pressure, PbTe(32j, 
SnTe[33], and GeTe|Ml- Provided that the Fermi sur- 
face is centered at time-reversal-invariant momenta, the 
Dirac-type relativistic k-p Hamiltonian gives a valid low- 
energy description of its band structure [26]. Moreover, 
if the pairing is odd under spatial inversion and fully 
gapped, the system is (almost) guaranteed to be a topo- 
logical superconductor according to our criterionpi 
In that case, our findings of unusual surface Andreev 
bound states apply directly. In addition, if the pairing 
symmetry of noncentrosymmetric superconductors (such 
as YPtBi[35]) has a dominant odd-parity component, 
their properties are often inherited from the centrosym- 
metric limit. We hope this work will encourage further 
explorations of unconventional surface Andreev bound 
states and stimulate the search for topological supercon- 
ductors. 

Note: Two recent studies [31 [TH] have calculated the 
surface spectral function numerically in Cu2:Bi2Se3 tight- 
binding models different from ours. The connection to 
our work on Majorana fermion surface Andreev bound 
state remains to be understood. 
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